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Abstract. We study the dynamics of a generalized inflationary model in which both the scalar field 
and its derivatives are coupled to the gravity. We consider a general form of the nonminimal derivative 
coupling in order to have a complete treatment of the model. By expanding the action up to the second 
order in perturbation, we study the spectrum of the primordial modes of the perturbations. Also, 
by expanding the action up to the third order and considering the three point correlation functions, 
the amplitude of the non-Gaussianity of the primordial perturbations is studied both in equilateral 
and orthogonal configurations. Finally, by adopting some sort of potentials, we compare the model 
at hand with the Planck 2015 released observational data and obtain some constraints on the model’s 
parameters space. As an important result, we show that the nonminimal couplings help to make 
models of chaotic inflation, that would otherwise be in tension with Planck data, in better agreement 
with the data. This model is consistent with observation at weak coupling limit. 
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1 Introduction 

The existence of an inflationary stage in the history of the early universe was proposed to solve some 
problems of the standard big bang cosmology, such as the flatness, horizon and relics problems. A 
successful inflationary paradigm can also provide a causal mechanism for generating the initial density 
perturbations needed to seed the formation of structures in the universe [1, 2, 3, 4, 5, 6, 7, 8]. A 
simple single field inflationary scenario (in which the early Universe was dominated by the potential 
energy of the scalar field dubbed inflaton) predicts the dominant mode of the primordial density 
fluctuations to be highly adiabatic, scale invariant and Gaussian [9]. But, recent observational data 
have detected a level of scale dependence in the primordial perturbations [10, 11]. On the other 
hand, although there is no direct signal for primordial non-Gaussianity in observation, Planck team 
has obtained some tight limits on primordial non-Gaussianity [12]. Some inflationary models also, 
predict a level of non-Gaussianity in the perturbation’s mode [9, 13, 14, 15, 16, 17, 18, 19, 20]. This is 
really an important and interesting point. Because, a large amount of information on the cosmological 
dynamics deriving the initial inflationary expansion of the Universe can be found in the primordial 
scalar non-Gaussianities [21, 22]. In this regard, the extended inflationary models, which can show 
the non-Gaussianity and scale dependence of the primordial perturbation, are more favorable. 

One successful class of the extended inflationary models is the one with a nonminimal coupling 
between the scalar field and the curvature scalar R. This coupling usually is shown by the lagrangian 
term as f{(f>)R. There are some motivations behind including an explicit nonminimal coupling in 
the action. It is necessary for the renormalization term arising in quantum field theory in curved 
space [23]. In considering the quantum corrections to the scalar field theory, this term would arises 
at the quantum level [24, 25]. With nonminimal coupling term, it is possible to have gravity as a 
spontaneous symmetry-breaking effect [26, 27]. NMG term allows to have an oscillating Universe [28]. 
Also, with this term, it is possible to solve the graceful exit problem of the old inflation model by 
slowing the false-vacuum expansion [29, 30]. So, it seems that the presence of nonminimal coupling 
between the scalar field and Ricci scalar is an useful extension of the theory and in this regard many 
authors have studied such models [17, 31, 32, 33, 34, 35, 36, 37, 38, 39]. 

It has been shown that the presence of NMC makes the inflation model more viable. For instance, 
the author of Ref. [40] has considered an inflation model with nonminimal coupling and non-canonical 
kinetic terms leading to an chaotic inflation model. It is shown that this chaotic inflation model 
is favored by the BIGEP2 observation. But this model is in tension with recent observational data 
released by Planck2015. The authors of Ref. [41] have considered an inflation model with a nonminimal 
coupling between the scalar field and gravity as {acj)^ + l3cj)^)R. They compared their model with 
Planck2015 data in Einstein frame to explore the viability of the model. It is shown that this model 
just for N = 65 and with some maximum values of is consistent with Planck2015 data. In 
Ref. [42] a nonminimally coupled inflaionary model with quadratic potential in Einstein frame has 
been considered. The authors of this paper have shown that their model for some values of the NMG 


- 1 - 




parameter is consistent with observational data. For other literature dealing with these issues we refer 
the reader to Refs. [43, 44, 45, 46, 47, 48, 49, 50, 51, 52]. 

Another extension of the inflation theory arises from considering a nonminimal coupling between 
the derivatives of the scalar field and the curvature, which leads to the interesting cosmological 
behaviors [53]. This coupling is usually given by the lagrangian term as V^cjiV'cj). 

Some authors have shown that during inflation, with a nonminimal coupling between the scalar field 
and gravity, the unitary bound of the theory violates [54, 55, 56]. However, it has been shown that the 
model with a nonminimal couplings between the kinetic term of the inflaton (derivatives of the scalar 
field) and the Einstein tensor preserves the unitary bound during inflation [57]. Also, the presence of 
nonminimal derivative coupling is a powerful tool to increase the friction of an inflaton rolling down its 
own potential [57]. Some authors have considered the model with this coupling term and have studied 
the early time accelerating expansion of the universe as well as the late time dynamics [58, 59, 60]. 
As an important extension of the nonminimal models, in addition to the standard coupling between 
the scalar field and curvature, it is possible also to have coupling between the derivatives of the scalar 
field and gravitational sector which makes the inflationary model more powerful to be a successful 
scenario. 

In this work, we consider a generalized inflationary model in which both the scalar held and its 
derivatives are coupled to the curvature. We also, consider an extension of the nonminimal derivative 
term in the lagrangian as F{(j))G{(j)). In this term, is a general function of (j) which for 
R ^(j), leads to the simple cases introduced earlier in literature. In section 2 we present some 
main equations of this generalized inhationary model. In section 3, we use the ADM metric and 
study the linear perturbations of the model. By expanding the action up to the second order in 
perturbation, and considering the 2-point correlation functions, we obtain the amplitude of the scalar 
perturbation and its spectral index. Also, by considering the tensor part of the perturbed metric, 
we obtain the tensor perturbation and its spectral index as well. In section 4, by expanding the 
action up to the cubic order in perturbation, we explore the non-linear perturbation in this model. 
By considering the 3-point correlation functions, we study the non-Gaussian modes of the primordial 
perturbations. In this section we obtain the amplitude of the non-Gaussianity in the equilateral and 
orthogonal conhguration and in ki = k 2 = ks limit. Then, in section 5, we test our generalized 
inhationary model in confrontation with the recently released observational data. To this end, we 
consider two types of nonminimal coupling parameter ( and several types of nonminimal 

derivative coupling function and also, several types of potential. By these functions we study the 
behavior of the perturbative parameters in the background of the Planck 2015 observational data. 
Our main goal is to see whether this non-minimal model has the potential to make models of chaotic 
inhation in better agreement with Planck 2015 data or not. 

In this regard, the model with non-minimal coupling between the scalar held and its derivative 
to gravity is a subset of Horndeski theory which satishes the necessary conditions for being a good 
theoretical inhationary model. Also, this model is observationally viable, since this model with several 
inhationary potential is consistent with Planck2015 observational data. As we can see from Planck 
team’s released paper, a simple inhationary model, with potential corresponding to (jP, (j>^ and (jy^ 
is not consistent with observation. In Ref. [52] it has been shown that a nonminimal model (a 
model with nonminimal coupling between the scalar held and gravity) is viable with potentials they 
have considered. However, their model has been studied in Einstein frame. We will show that 
adding a nonminimal coupling between the gravity and derivatives of the scalar held leads to a viable 
inhationary model even in Jordan frame. 

The main properties of the cosmological perturbations are described by the tensor-to-scalar ratio 
and the scalar spectral index. In this regard, several observational teams are trying to hnd some con¬ 
straints on these perturbative parameters. WMAP team has found the constraints r < 0.13 and Us = 
0.9636 ± 0.0084 from the combined WMAP9-|-eGMB-|-BAO-|-HQ data [61]. The constraints obtained 
by the joint Planck 2013-1-WMAP9-I-BAO data are as r < 0.12 and = 0.9643 ± 0.0059 [62]. Then, 
in 2014, BIGEP team found a surprising result. The joint Planck 2013-1-WMAP-I-BIGEP2-I-BAO 
data has set the constraints r = 0.2096j'jQ Qggg and Us = 0.9653 ± 0.0129 on the perturbative pa¬ 
rameters [63, 64]. This large value of the tensor-to-scalar ratio, which shows the large amplitude of 
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the gravitational wave modes generated during inflation, was a really challenging result. But Planck 
collaboration has performed a new study on the polarized dust emission in our galaxy and has shown 
that the part of the sky observed by BICEP2 is contaminated significantly by galactic-dust emission. 
So, the polarized dust emission might contribute a significant part of the B-mode detected by BI- 
CEP2. The Planck collaboration released the constraints r < 0.099 and 0.9652 ± 0.0047 from Planck 
TT, TE, EE-flowP -fWP data (note that Planck TT, TE, EE-l-lowP refer to the combination of the 
likelihood at / > 30 using TT, TE, and EE spectra and the low I multipole temperature polarization 
likelihood) [10, 11, 12]. In this regard, to compare our model with observational data, we study the 
behavior of the tensor-to-scalar ratio versus the scalar spectral index in the background of the Planck 
TT, TE, EE-flowP data and deduce some constraints on the parameters space of the model. We 
also analyze the non-Gaussianity feature of the model by studying the behavior of the orthogonal 
configuration versus the equilateral configuration in the background of the observational data. 


2 The Setup 


The four-dimensional action for a cosmological model where a scalar field is nonminimally coupled 
to the Ricci scalar, in the presence of a nonminimal derivative coupling between the scalar field and 
gravity, is given by the following expression 


5 = 


d'^xyf—g 




( 2 . 1 ) 


where, R is the Ricci scalar, ()) is a scalar field (inflaton) with the potential V{(j)) and f{(j)) and J^(<(') 
are general functions of the scalar field. Also, in the nonminimal derivative term of the action (2.1), 
= R^u — ^g^yR is the Einstein tensor. Note that, in some papers there is a coefficient ^ in front 
of the nonminimal derivative term, where the constant parameter k* has dimension of length-squared. 
In this paper, we absorb such a coefficient in function By assuming the ERW metric, variation 

of the action (2.1) with respect to the metric leads to the Eriedmann equation of the model as 

(1 + 2K^f) = y ^ ^ (2.2) 

where, a dot refers to a time derivative of the parameter and a prime denotes a derivative with respect 
to the scalar field. By varying the action (2.1) with respect to the scalar field we get the following 
equation of motion 

- 1 -t- + {l2R'HH + IST'H^ - Si?) ^ + 6fR - W = 0. (2.3) 

Now, by defining the slow-roll parameters as e = — and ry = "we find these parameters 

in our setup as follows 

e =- - -(2.4) 

1 -f 2^2/ + 


^ A A 2K^f'^ -f 

^ i7e(l -f 2^2/ -p Kjr'<^ 2 ) ITt (1 -P 2k 2/ -P Kjr'(^ 2)2 ’ 

where parameter A is defined as 

A = —At - — - 3k^H-H- — + H- — ■ 

2iJ2 H ^ H H<l)\ H H ) 


(2.5) 


( 2 . 6 ) 


During the inflationary era, the evolution of the Hubble parameter is so slow, so that in this era the 
conditions e <C 1 and rj 1 are satisfied. As one of these two slow-varying parameters reaches unity, 
the inflation phase terminates. 













The number of e-folds during inflation which is defined as 


ftf 

N= Hdt, (2.7) 

^hc 


in our setup and within the slow-roll limit {(p |3i7^| and (j)^ <C V{(j))) takes the following form 

r<t>f - l) 


N ~ 


' 4*hc 


V' - 6f'R 


-dp, 


( 2 . 8 ) 


where ptc shows the value of the inflaton at the horizon crossing of the universe scale and p / denotes 
the value of the field when the universe exits the inflationary phase. In the next section, we study 
the linear perturbation of the model which is the key test of any inflationary model. In this regard, 
we calculate the spectrum of perturbations produced due to quantum fluctuations of the fields about 
their homogeneous background values. 


3 Linear Perturbation 


In this section, we study the linear perturbation arising from the quantum behavior of both the scalar 
field p and the space-time metric, around the homogeneous background solution. To this end, 
we should expand the action (2.1) up to the second order in small fluctuations. In this regard, it is 
convenient to work in ADM metric formalism [65] 

= —N'^dt^ -\- hij [dx^ -\- N^dt) [dx^ -\- dt ), (3-1) 

where N is the lapse function and N'^ is the shift vector. To obtain a general perturbed form of 
the metric (3.1), we should expand the shift and laps functions as N'^ = B'^ = 6'^^djB -|- u® and 
N = 1 -I- respectively. $ and B are 3-scalar and u® is a vector satisfying the condition = 
0 [66, 67]. Note that it is sufficient to compute N or up to the first order. This is because the 
second order perturbation is multiplied by a factor which is vanishing using the first order solution. 
Since the third order term is multiplied by a constraint equation at the zeroth order obeying the 
equations of motion, the contribution of this term also vanishes [9, 15, 68]. hy should be written as 
hij = a? [(1 — 2'i>)Sij 2Tij], where 'k is the spatial curvature perturbation and Tij is a spatial shear 
3-tensor which is symmetric and traceless. So, the perturbed metric (3.1) takes the following form 

= _(i + 2^)dt‘^ + 2a(t)B, dtdx^ + a^{t) [(1 - 2^-)^^ -b 27^^]. (3.2) 


In order to study the scalar perturbation of the theory, we work within the uniform-field gauge with 
5p = 0 and also the gauge Tij = 0. So, by considering the scalar part of the perturbations at the 
linear level, we rewrite the perturbed metric (3.2) as [66, 67, 69] 

ds^ = —(1 -b 2$)dt^ -b 2a{t)Bj dtdx'' -b aP'{t){l — 2^)5ij. (3.3) 

Replacing the metric (3.3) in action (2.1) and expanding the action up to the second order in pertur¬ 
bations gives 


82= i dtd^xa^ 


-3(k-" -bZ-brJ")^ - 


2(k " + / + , 1 


+ -L(2{n-^ + f + P^B')i> 


-{2k-'^H{ 1 + K^f) + 2pf + 6p'^T')<i>'jd^B + 3(2 k-'^H{ 1 + n^f) + 2pf + Qp^T''^ 
- (3K"2ij2(l -b K^f) - ^p^ + QHpf + -b 


K” + / - 


( a «')2 


(3.4) 
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In this gauge, the equations of motion for $ and B, obtained by varying the action (3.4), lead 
to the following constraint equations 


$ = 2 - 


K ^ + / + (j?P 


+ K^f) + 20/' + 




(3.5) 


d^B 


= 34'- 


+ f + 


+ K^f) + 20/' + 602 J-') 

2(9«;-2ij2(i + pf) _ 1^2 iSH^f + 5417202jr/) 
3(2K-2i7(l + k 2/) + 20/' + 602j-/) 


-dH 


$. 


By substituting equation (3.5) in action (3.4) and integrating by parts, we reach 


/ 


S '2 = / dtPxPW 


4-2 _ ^(a^)2 


where the parameters W and (which is known as the sound velocity) are defined as 

, 2 


W = -4 


(9^-2172(1 + Pf ) - |02 + I8i70/' + 5471202 

(2K-27r(l + k2 /) + 20/' + 602 

+3 ^ + f + (/>^p'j , 


(3.6) 


(3.7) 


(3.8) 


c? = 


?(2K-2i7(l + Pf) + 20/' + 60^J-') + / + (p^pyH - (2 k-'^H{ 1 + Pf) + 20/' + 

(^-2 + f_ ^ 4(2k-^H{ 1 + Pf) + 20/' + 602 J-') (-^-2 + / + 02 J-') ^- - - 1 

-2(.-2 + / + + I f ^ ^ 


3 > 


1 -1 


+602_ 4j'^-2 ^2 ( 9 ^- 2772(1 + pf) - ^02 18770/' + MH^P^P 


(3,9) 


respectively. To see more details of driving this type of equations, one can refer to Refs. [15, 16, 
22, 70]. 

Now we calculate the quantum perturbations of 4^. To this end, by variation of the action (3.7), 
we find the equation of motion of 4^ as 


^+13i7+J^|4' + c2^^ = 0. 


(3.10) 


The solution of the above equation, up to the lowest order of the slow-roll approximation, is given by 
the following expression 

jfXp-iclkT 

^ = 3 (1 + tc2fcr) . (3.11) 

2 c| y/PW 
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To study the power spectrum of the curvature perturbation, we should compute the two point corre¬ 
lation function in our setup. We find the two-point correlation function, sometimes after the horizon 
crossing, by obtaining the vacuum expectation value of dl at r = 0 (corresponding to the end of 
inflation) 

27r2 

(0|«'(0, ki)^(0, k2)|0) = (27r)353(ki -b k2)-p-A, (3.12) 

where As is called the power spectrum and is given by 


H'^ 


(3.13) 


The scalar spectral index of the perturbations at the Hubble crossing is defined as 


ns - 1 


d In As 


din fc 


Cs k—aH 


which in our setup takes the following form 


Us — I = —2e — 


5k^/'0 
2H{1 + K^f) 


1 dln(e-b ^ 1 d\ncs 

H dt H dt 


(3.14) 


(3.15) 


Any deviation of rig from the unity shows that the evolution of perturbations is scale dependent. 

Other important parameters in an inflationary model are amplitude of the tensor perturbation 
and its spectral index. To study the tensor perturbation we use the tensor part of the perturbed 
metric (3.2). We can write 7ij, in terms of the two polarization tensors, as follows 

7-,=r+e++rxer^, (3.16) 

where e]J’ ' are symmetric, traceless and transverse. The normalization condition imposes the con¬ 
straints 

e?’"^(k)e(+’")(-k)* = 2, (3.17) 

e«(k)e(;)(-k)*=0, (3.18) 

and the reality condition gives 

e(+’^)(-k) = (e(;’")(k))*. (3.19) 

Now, we can write the second order action for the tensor mode (gravitational waves) as follows 


St 


dt d^x a^Wr 


fl - ^{dr+f + fi - 


(3.20) 


where Wt and Cy are defined with the following expressions respectively 




(3.21) 


C'-p — 


k- 2 -b / -b 4>'^r 


(3.22) 


By following the strategy as applied for the scalar mode case, we can find the amplitude of the tensor 
perturbations as follows 
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which, by using the definition of the tensor spectral index 


TIT = 

leads to the following expression for tit 

ttt = —2e — 


dlnAx 

dink ’ 




2i7(l + /t2/)' 


(3.24) 


(3.25) 


The ratio between the amplitudes of the tensor and scalar perturbations (tensor-to-scalar ratio) is 
another important inflationary parameter which in our setup is given by 


^ = rll = i6c (cl 1 

At H 4i?(l + K^f) ) ■ 


(3.26) 


4 Nonlinear Perturbations and Non-Gaussianity 

Non-Gaussianity of the primordial density perturbations is another important aspect of an infla¬ 
tionary model. To explore the non-Gaussianity of the density perturbations, one has to study the 
nonlinear perturbation theory. The two-point correlation function of the scalar perturbations gives 
no information about the non-Gaussianity of the model. The non-Gaussian property shows itself up 
in the three-point correlation function. This is because for a Gaussian perturbation, all odd n-point 
correlators vanish and the higher even n-point correlation functions can be expressed in terms of sum 
of products of the two-point functions. To calculate the three-point correlation function, we should 
expand the action (2.1) up to the cubic order in the small fluctuations around the homogeneous back¬ 
ground solution. The obtained cubic terms, lead to a change both in the ground state of the quantum 
field and non-linearities in the evolution. By expanding the action (2.1) up to the third order in per¬ 
turbation, we obtain a complicated expression which can be seen in the Appendix 6. We use equation 
(3.5) to eliminate the perturbation parameter <i> in the expanded action. Then by introducing the 
new parameter X as 


B = 


2(K-" + / + r^')^ 


-k 


a^X 


and 


2k-'^H{ 1 + K^f) + 2^f' + 6(^2jr/ ^-2 j ^2jr/ ’ 


d^X = W4, 


(4.1) 

(4.2) 


we can rewrite the third order action, up to the leading order, as follows 



1 + «V 




(4.3) 


Now, we are in the position to calculate the three point correlation function by using the inter¬ 
acting picture. In this picture, we obtain the vacuum expectation value of the curvature perturbation 
41 for the three-point operator in the conformal time interval between the beginning and end of the 
inflation as [9, 22, 70] 

(4'(ki) 4'(k2) ^(ks)) =-z J (iTa(0|[4'(T/,ki)4'(T/,k2)4'(r/,k3),i7i„t(T)]|0), (4.4) 
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where the interacting Hamiltonian, Hint, is equal to the lagrangian of the third order action. To solve 
the integral of equation (4.4), we can approximate the coefficients in the brackets of the lagrangian 
(4.3) to be constants since these coefficients would vary slower than the scale factor. By solving the 
integral, we find the three-point correlation function of the curvature perturbation in the Fourier 
space as 


where 


where 


and 


^^(ki, k2, ka) — 


Vi + T 


2 

'k 


i>j 

^ K 

S. = 


3^3(ki+k2-t-k3)H^(ki,k2,k3), 

(4.5) 

(27r)'‘A^ 

nL^.» *( '■ “”)■ 

(4.6) 

). Also, the parameter is defined as 



(4.7) 

2 ^2_L V ^3 

(4.8) 

i>j iiij 

(4.9) 

(ki k2 ksf 

(4.10) 


(4.11) 


As it can be seen from equation (4.6), a three-point correlator depends on the three momenta ki , k 2 
and ks- To satisfy the translation invariance, these momenta should form a closed triangle, meaning 
that the sum of these momenta should be zero (fci + k 2 + k^ = 0). Also, satisfying the rotational 
invariance makes the shape of the triangle important [21, 71, 72, 73, 74]. Depending on the values 
of momenta, there are several shapes and each shape has a maximal signal in a special configuration 
of triangle. One of these shapes is a local shape [75, 76, 77, 78] which has a maximal signal in the 
squeezed limit (fca <C fci — ^ 2 ). Another shape which is corresponding to the equilateral triangle [79], 
has a signal at fci = /c 2 = fca. A linear combination of the equilateral and orthogonal templates, which 
is orthogonal [80] to equilateral one, gives a shape corresponding to folded triangle [81] with a pick in 
ki = 2k2 = 2^3 limit. We mention that, the orthogonal configuration has a signal with a positive peak 
at the equilateral configuration and a negative peak at the folded configuration. In order to measure 
the amplitude of the non-Gaussianity we define the dimensionless parameter called “nonlinearity 
parameter”, as follows 


^ 10 


(4.12) 


In this paper, we study the non-Gaussianity in the equilateral and orthogonal configurations. To this 
end, we should find in these configurations. In this regard, we should estimate the correlation 
between two different shapes. So, following Refs. [82, 83, 84], we define the following quantity 






(4.13) 


where, ^ and 


^■(H^H^) = / dkidk2dkzB,i,{ki,k2,k3)B't^{ki,k2,k^)w 


( 4 . 14 ) 







with w = ■ The region of integration is 0 < fci < oo, 0 < |^ < 1 and < 1- 

We note that two shapes satisfying the condition | €{3^3'^) |~ 0, are almost orthogonal. Now, we 
introduce a shape 5®?“® as [20, 83] 

=-^( 351 - 52 ). (4.15) 

By using this equation and equation (4.14) we can show that the foilwing shape is orthogonal to (4.15) 

- '^2) + 35i - 52) , (4.16) 

where /3 ~ 1.1967996. Now, the bispectrum (4.7), can be expressed in terms of the equilateral and 
orthogonal basis as 

Sq, = Cl + C 2 , (4.17) 

where the parameters Ci and C 2 are given by the following expressions 


Cl 


13 

12 



2 + 3/3 


(4.18) 


and 


C 2 


14 - 13/3 
12 



(4.19) 


By using equations (4.12)-(4.19), we can obtain the amplitude of the non-Gaussianity in the equilateral 
and orthogonal configuration as follows 


requi _ 

^ NL 


130 

36Etifc? 



3/3 


^equi 


(4.20) 


and 


northo 
^ NL 


140 - 130^ 
36 



(4.21) 


Since the equilateral shape has a maximal signal at ki = k 2 = k^ limit and also, an orthogonal 
configuration has a positive peak at this limit, so we rewrite equations (4.20) and (4.21) in this limit 
as 


requi _ 

J NL 


325 


1 / 1 


24 Vc? 


-12 + 3/3 


(4.22) 


and 


Jnl 9 V 4 ^ 6/ 8 V 


1 f 1 


(4.23) 


So far, we have obtained the main equations of our setup. In the following section, we examine 
our model by using the recently released observational data by Planck Collaboration and find some 
constraints on the model’s parameters space. 


5 Observational Constraints 

In this section we are going to find some observational constraints on the parameters space of the 
model in hand. To this end, first we define the extension of the nonminimal derivative term as 
~ 2 n^”’ which for n = 1 leads to the usual nonminimal derivative case (as mentioned in the 
Introduction). We set the values I, 2, 3 and 4 for n. Another general function in this setup is 
/((/)), which shows the nonminimal coupling between the scalar field and the curvature scalar. For 


- 9 - 



















this function we consider two cases: f(^) ~ and /(</>) ~ . The last general function is the 

potential term in the action. We adopt various types of potential to explore: the linear potential 
as V{(j)) ~ (j) [85], the quadratic potential as V{(j)) ~ , the cubic potential as V{(j)) ~ , the 

quartic potential as V{(1)) ~ 4>‘^, the axion monodromy’s motivated potential as V{(j)) ~ ^3 [86] and 
the exponential potential as ^ With this specification of the general functions, we can 

analyze the model numerically and find some constraints on the nonminimal coupling parameter 
For this purpose, we solve the integral of equation (2.8) by adopting the mentioned functions of 
the scalar field. By solving this equation, we find the value of the inflaton at the horizon crossing 
in terms of the e-folds number. By substituting (j)hc in equations (3.15), (3.26), (4.22) and (4.23), 
we can find the scalar spectral index, tensor-to-scalar ratio, the amplitude of the orthogonal and 
equilateral configuration of the non-Gaussianity in terms of the e-folds number. Then we can analyze 
these parameters numerically to see the observational viability of this setup in confrontation with 
recently released data. In plotting the figures and analyzing the model we have re-scaled all constant 
parameter to unity with the choice G = c = 1. In this regard, we study the behavior of the tensor-to- 
scalar ratio versus the scalar spectral index and the orthogonal configuration versus the equilateral 
configuration. The observational parameters are defined at fcp = 0.002Mpc“^ where subscript 0 refers 
to the value of k when it left the Hubble radius during inflation. The results are shown in figures 
1-8. We have found that an inflationary model with a nonminimal coupling between the scalar filed 
and Ricci scalar and also a nonminimal derivative coupling, in some ranges of nonminimal coupling 
parameter ^ is consistent with Planck20I5 dataset. In fact, these nonminimal couplings help to make 
models of chaotic inflation, that would otherwise be in tension with Planck data, in better agreement 
with the data. The ranges of nonminimal coupling are shown in tables I-IV. Note that, some other 
authors have considered inflationary model with nonminimal coupling term and studied their model 
numerically. We can compare the result of our setting with their results here. For example, it is 
shown in Ref. [40] that the presence of nonminimal coupling and non-canonical kinetic terms lead to 
an chaotic inflation model, which is favored by the BICEP2 observation but is in tension with recent 
observational data. Here, we have shown that the presence of nonminimal coupling between the scalar 
field and Ricci scalar and also a general nonminimal derivative coupling between the scalar field and 
gravity make the model observationally more viable. This model with these couplings is consistent 
with Planck20I5 TT, TE, EE-|-lowP data. The authors of Ref. [41] have considered a nonminimal 
coupling between a scalar field and gravity as (ai/)^ -I- in Jordan frame. They perform a frame 

transformation to Einstein frame and compare their model with Planck2015 data. Their model just 
for TV = 65 and with some maximum values of Us is consistent with observational data and for 
N = 50 and 60 is observationally disfavored. In our setup, we have considered nonminimal coupling 
function (j)^ and separately. Our numerical analysis shows that the presence of NMC and NMDC 
makes the model observationally viable for both TV = 50 and = 60 in the Jordan frame (note that 
although we don’t present here the results with N = 50, but we have analyzed the model with this 
value of number of e-folds and have found the viability of the model). In Ref. [42] a nonminimally 
coupled inflaionary model with quadratic potential in Einstein frame is considered in the background 
of Planck20I5 observational data. Their numerical analysis shows that this model for some value of 
NMC parameter gives Ug ~ 0.9652. Interestingly our setup also, in some values of ^ gives Ug ~ 0.9652. 
Nevertheless, we should emphasize that for some other choices for model parameters, this would be 
not the case necessarily. 

Our numerical analysis shows that this nonminimal model with small values of the NMC param¬ 
eter is consistent with observational data and so this model in Jordan frame in weak coupling limit has 
cosmological viability. In week coupling limit our setup is applicable to several inflationary potentials. 
In this limit and in the ranges obtained by numerical analysis, the scalar spectral index is red tilted. 
Another point is that in our setup, in some ranges of NMC parameter, it is possible to have large 
non-Gaussianity, as can be seen in figures. Note that, to obtain the constraints on we have focused 
on those ranges of the NMC parameter that lead to the values of the observable quantities allowed 
by the Planck 2015 observational data. For instance, we have focused on the values of ^ leading to 
0.9522 < Ug < 0.977 allowed by Planck20I5 TT, TE, EE-flowP data. In the same way, we focused on 
the values of the NMC parameter leading to —147 < < 143 allowed by Planck20I5 TTT, EEE, 
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Figure 1. Tensor-to-scalar ratio versus the scalar spectral index for an inflationary model in which both the 
inflaton and its derivative are nonminiinally coupled to gravity, in the background of the Planck2015 TT, 
TE, EE+lowP data (left panel) and the amplitude of the orthogonal configuration of the non-Gaussianity 
versus the amplitude of the equilateral configuration for an inflationary model in which both inflaton and its 
derivative are nonminimally coupled to gravity, in the background of Planck2015 TTT, EEE, TTE and EET 
data (right panel). Figure is plotted with n = 1, f{(j}) ~ S’Hd various types of the potential, for = 60 


Table 1. The ranges of ^ in which the values of the inflationary parameters r and Ua and also, f^L^° and 
with n = 1, are compatible with the 95% CL of the Planck2015 dataset. 


V 

r - ns 

r - Tis 

fortho fG-qm 

JNT. ~ Jmt. 

fOrtho fGqm 

JNT. - Jmt. 


fW ~ 

fW ~ 

f{<f>) ~ 



€ < 0.091 

C < 0.0912 

C < 0.090 

0.008 < 5 < 0.091 


€ < 0.088 

? < 0.09 

0.01 < € < 0.089 

0.015 < € < 0.0895 


0.012 < € < 0.09 

€ < 0.08 

0.011 < c < 0.093 

0.021 < C < 0.082 


0.014 < € < 0.10 

0.001 < € < 0.093 

0.013 < € < 0.10 

0.01 < C < 0.098 

2 

(ji 3 

€ < 0.081 

€ < 0.082 

0.015 < C < 0.073 

0.01 < C < 0.109 


0.016 < C < 0.108 

0.078 < C < 0.080 

0.094 < C < 0.103 

0.017 < C < 0.097 


TTE and EET data. Also, equations (3.15)), (3.26), (4.22) and (4.23) show that the scalar spectral 
index, tensor-to-scalar ratio and the amplitudes of non-Gaussianity have nonlinear dependence on 
the NMC parameter The function /((/>) has an explicit non-miniinal coupling dependence. In the 
mentioned equations, there are different powers of the function /(</>) and therefore different powers of 
the nonminimal coupling parameter. Because of complexity of equations we have avoided to rewrite 
the equations in terms of the NMC parameter However, we note that the effects of these nonlinear¬ 
ities in ^ have been included in our analysis and graphs. Since quantities such as the scalar spectral 
index, tensor-to-scalar ratio and the amplitudes of non-Gaussianity have nonlinear dependence on the 
NMC parameter any small change in the values of ^ has considerable effect on the values of these 
quantities. In fact, there is a nonlinear dependence on model parameters that really demands four 
figures of accuracy in quantities such as ^ to get two figures of accuracy in observable quantities such 
as rig. So, we have tried to be more accurate and found the values of ^ with higher precision. 
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Figure 2. Tensor-to-scalar ratio versus the scalar spectral index for an inflationary model in which both 
inflaton and its derivative are nonminimally coupled to gravity, in the background of the Planck2015 TT, TE, 
EE+lowP data (left panel) and the amplitude of the orthogonal configuration of non-Gaussianity versus the 
equilateral configuration for an inflationary model in which both inflaton and its derivative are nonminimally 
coupled to gravity, in the background of Planck2015 TTT, EEE, TTE and EET data, (right panel). Figure 
is plotted with n = 1, /(<(>) ~ and various types of potential, for N = 60. 



Figure 3. Tensor-to-scalar ratio versus the scalar spectral index for an inflationary model in which both 
inflaton and its derivative are nonminimally coupled to gravity, in the background of Planck2015 TT, TE, 
EE-|-lowP data (left panel) and the amplitude of the orthogonal configuration versus the equilateral configu¬ 
ration of non-Gaussianity for an inflationary model in which both inflaton and its derivative are nonminimally 
coupled to gravity, in the background of Planck2015 TTT, EEE, TTE and EET data (right panel). Figure is 
plotted with n = 2, ~ and various types of potential, for N = 60. 

6 Summary 

In this paper, we have studied the dynamics of a generalized inflationary model in which both the 
scalar field and its derivatives are nonminimally coupled to gravity. In order to have a complete 
treatment of the model, we have considered an extension of the nonminimal derivative coupling in 
the lagrangian as F{(j))G{(j)) where is a general function of <j). We have firstly obtained 
the main equations of this generalized inflationary model. Then, using the ADM formulation of the 
metric, we have studied the linear perturbations of the model. We have expanded the action up 
to the second order in perturbation and by using the 2-point correlation function, we have derived 
the amplitude of the scalar perturbation and its spectral index in our setup. To study the tensor 
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Table 2. The ranges of ^ in which the values of the inflationary parameters r and ris and also, f^L^° and 
fNL^y with n = 2, are compatible with the 95% CL of the Planck2015 dataset. 


V 

r - ns 

r - Us 

fortho 


fortho fSquz 

JNT. ’ J NT. 


fW ~ 

fW ~ 

n*) ~ 




0.021 < € < 0.090 

0.019 < 5 < 0.09 

0.020 < € 

< 0.088 

0.021 < € < 0.097 


€ < 0.109 

0.01 < 5 < 0.055,0.06 < € < 0.114 

o 

b 

A 

< 0.89 

0.014 < 5 < 0.109 


0.022 < € < 0.081 

0.009 < € < 0.086 

V 

00 

o 

b 

< 0.098 

€ < 0.086 


€ < 0.098 

0.010 < € < 0.09 

V 

o 

b 

< 0.107 

0.0108 < € < 0.098 

2 
(j) ^ 

0.021 < € < 0.076 

0.018 < € < 0.084 

V 

o 

CN 

o 

b 

< 0.086 

0.011 < € < 0.084 

e-^ 

€ < 0.086 

0.01 < C < 0.081 

0.028 < € 

< 0.089 

0.01 < ^ < 0.104 



^qui 

■'NL 



-V~<|) ,-v~l|l^ ,-V~l|l’ , 


-V~<|) ,-v~l|l^ ,-V~l|l’ , 

' V~l|)''{2/3} ,-V~e 


' V~l|)''{2/3} ,-V~e 


Figure 4. Tensor-to-scalar ratio versus the scalar spectral index for an inflationary model in which both 
inflaton and its derivative are nonminimally coupled to gravity, in the background of Planck2015 TT, TE, 
EE+lowP data (left panel) and tThe amplitude of the orthogonal configuration versus the amplitude of the 
equilateral configuration of non-Gaussianity for an inflationary model in which both inflaton and its derivative 
are nonminimally coupled to gravity, in the background of Planck2015 TTT, EEE, TTE and EET data (right 
panel). Figure is plotted with n — 2, f(4>) ~ and various types of potential, for N = 60. 


Table 3. The ranges of ^ in which the values of the inflationary parameters r and ris and also, ° and 
with n = 3, are compatible with the 95% CL of the Planck2015 dataset. 



r - ns 

r - 2 

T'S 


fOrtho fBqm 

J NT. ' J NT. 

fortho 

J NT. 

fCqui 

^NL 


fW ~ 

fW ~ 


/(<A) ~ e<A^ 

fW ~ 



0.020 < C < 0.089 

0.017 < € 

< 

0.089 

0.018 < C < 0.107 

0.011 

< e 

< 0.089 


0.015 < C < 0.084 

0.012 < C 

< 

0.090 

0.011 < C < 0.098 

0.012 

< C 

< 0.090 

<i>^ 

0.018 < C < 0.085 

V 

CD 

O 

O 

b 

< 

0.083 

0.02 < C < 0.079 

0.019 

< e 

< 0.080 


0.007 < C < 0.108 

0.032 < C 

< 

0.105 

0.029 < C < 0.081 

0.024 

< C 

< 0.871 

2 

<j> 3 

0.04831 < C 

0.011 < c 

< 

0.082 

0.086 < € < 0.096 

0.007 

< C 

< 0.095 


e-* 


0.013 < S < 0.088 0.043 < 5 < 0.078 0.025 < S < 0.077 0.021 < 5 < 0.081 
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Figure 5. Tensor-to-scalar ratio versus the scalar spectral index for an inflationary model in which both 
inflaton and its derivative are nonminimally coupled to gravity, in the background of Planck2015 TT, TE, 
EE+lowP data (left panel) and the amplitude of the orthogonal configuration versus the amplitude of the 
equilateral configuration of non-Gaussianity for an inflationary model in which both inflaton and its derivative 
are nonminimally coupled to gravity, in the background of Planck2015 TTT, EEE, TTE and EET data (right 
panel). Figure is plotted with n = 3, f{4>) ^ and various types of potential, for N = 60. 



Figure 6. Tensor-to-scalar ratio versus the scalar spectral index for an inflationary model in which both 
inflaton and its derivative are nonminimally coupled to gravity, in the background of Planck2015 TT, TE, 
EE-|-lowP data (left panel) and the amplitude of the orthogonal configuration versus the amplitude of the 
equilateral configuration of non-Gaussianity for an inflationary model in which both inflaton and its derivative 
are nonminimally coupled to gravity, in the background of Planck2015 TTT, EEE, TTE and EET data (right 
panel). Figure is plotted with n = 3, /((?i) ~ and various types of potential, for N = 60. 


perturbation we have used the tensor part of the perturbed metric and we have calculated the tensor 
perturbation and its spectral index too. The ratio between the amplitude of the tensor and scalar 
perturbations has been obtained in our setup. After that, in order to seek the non-Gaussian feature 
of the primordial perturbations, we have studied the non-linear theory. To study the non-linear 
perturbation, one has to calculate the three point correlation functions. In this regard, we have 
expanded the action of the model up to the third order in perturbation. Then, by using the interacting 
picture we have computed the three point correlation function and the nonlinearity parameter in our 
model. Also, by introducing the shape functions as 5®^"* and we have found the amplitude of 

the non-Gaussianity in the equilateral and orthogonal configurations. Since, both the equilateral and 
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Figure 7. Tensor-to-scalar ratio versus the scalar spectral index for an inflationary model in which both 
inflaton and its derivative are nonminimally coupled to gravity, in the background of Planck2015 TT, TE, 
EE+lowP data (the left panel) and the amplitude of the orthogonal configuration of non-Gaussianity versus the 
amplitude of the equilateral configuration for an inflationary model in which both inflaton and its derivative 
are nonminimally coupled to gravity, in the background of Planck2015 TTT, EEE, TTE and EET data (right 
panel). Figure is plotted with n = 4, f{4>) ^ and various types of potential, for N = 60. 



Figure 8. Tensor-to-scalar ratio versus the scalar spectral index for an inflationary model in which both 
inflaton and its derivative are nonminimally coupled to the gravity, in the background of Planck2015 TT, TE, 
EE-|-lowP data (left panel) and the amplitude of the orthogonal configuration of non-Gaussianity versus the 
amplitude of the equilateral configuration for an inflationary model in which both inflaton and its derivative 
are nonminimally coupled to the gravity, in the background of Planck2015 TTT, EEE, TTE and EET data 
(right panel). Figure is plotted with n = 4, ^ (,4>* and various types of potential, for N = 60. 


orthogonal configuration have peek in the limit ki = k 2 = ^ 3 , we have focused in this limit. After 
finding the important perturbation parameters, we have analyzed our model with recently released 
Planck 2015 observational data. To this end, we had to specify general functions of the scalar field in 
the model. For the the extension of the nonminimal derivative term, we have adopted the function 
to be as ~ with n = 1,2, 3,4. We have also defined two functions for the nonminimal 

coupling between the scalar field and the curvature scalar as f ^ and / ~ Then, by adopting 
several types of potential as linear, quadratic, cubic, quartic potential, a potential corresponding to 
03 and an exponential potential, we have analyzed this inflationary models numerically and compared 
our model with the observational data of Planck2015. In this regard, we have studied the behavior of 
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Table 4. The ranges of ^ in which the values of the inflationary parameters r and ris and also, f^L^° and 
/nl'’ with n = 4, are compatible with the 95% CL of the Planck2015 dataset. 


V 

r ~ Tig 

r - ns 

fortho fBqui 

J NT. ~ J NT. 

fOrtho fBqui 

J NT. ' J NT. 




/(0) ~ €0^ 

/(0) ~ €0^ 

<!> 

C < 0.084 

C < 0.021,0.025 < € < 0.084 

C < 0.113 

0.016 < C < 0.093 

<1,^ 

C < 0.083 

C < 0.040,0.059 < C < 0.112 

0.125 < C < 0.090 

0.012 < C < 0.092 


C < 0.082 

0.019 < 5 < 0.087 

C < 0.106 

0.007 < C < 0.109 


0.075 < ^ < 0.105 

0.019 < C < 0.081 

0.036 < C < 0.079 

0.019 < C < 0.0.935 

2 
(f) ^ 

0.043 < € < 0.112 

€ < 0.086 

0.009 < 5 < 0.099 

0.008 < C < 0.098 


0.012 < € < 0.099 

€ < 0.109 

0.004 < € < 0.104 

0.013 < € < 0.096 


the tensor-to-scalar ratio versus the scalar spectral index in the background of the Planck2015 TT, 
TE, EE+lowP data and obtained some constraints on the nonminimal coupling parameter. We have 
also analyzed the non-Gaussinaty feature of the primordial perturbations numerically by studying 
the behavior of the orthogonal configuration versus the equilateral configuration in the background of 
the Planck2015 TTT, EEE, TTE and EET data. In this paper, we have shown that this generalized 
inflationary model is observationally viable in some ranges of its parameters space. As an important 
result, we have shown that the nonminimal couplings help to make models of chaotic inflation, that 
would otherwise be in tension with Planck data, in better agreement with the data. 

In summary, as an important result, our numerical analysis shows that this non-minimal model, 
in some ranges of the non-minimal coupling parameter and with small values of this parameter (the 
weak coupling limit), is consistent with observation. Therefore, this model in weak coupling limit has 
cosmological viability. Also, in the weak coupling limit this setup is applicable to several inflationary 
potentials. In this limit and in the ranges obtained by numerical analysis, the scalar spectral index is 
red tilted. Another point is that in our setup, in some range of the NMC parameter, it is possible to 
have large non-Gaussianity. 
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Appendix: The expanded action up to third order 


S3 = I dtd^xa^ 
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